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We discuss the effects of imperfect photon detectors suffering from loss and 
noise on the reliabihty of linear optical quantum computers. We show that for 
a given detector efficiency, there is a maximum achievable success probability, 
and that increasing the number of ancillary photons and detectors used for 
one controlled sign flip gate beyond a critical point will decrease the prob- 
ability that the computer will function correctly. We have also performed 
simulations of some small logic gates and estimate the efficiency and noise 
levels required for the linear optical quantum computer to function properly. 
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Theoretically quantum computers can perform some types of calculations much faster 
than classical computers but the technological difficulties of manipulating quantum 
information have so far prevented researchers from constructing a quantum computer which 
is able to perform useful tasks. One of the first proposals for the construction of a quantum 
computer envisioned a single photon traveling through a network of beam splitters and phase 
shifters. Unfortunately for this single photon computer to encode and manipulate n qubits 
requires on the order of 2" optical modes, beam splitters, phase shifter, and photon detectors. 
The durability of the optical qubits and the simplicity of the linear optical components 
make this an attractive design for building small (a few qubits) quantum computers, but 
the exponential growth in the number of components makes this sort large scale quantum 
computer economically impractical. 

This problem of exponential growth seems to have been eliminated by Knill, Lafiamme 
and Milburn in [0-^]. In their proposal (outlined below) a single qubit is represented by 
the presence of a photon in one of two optical modes. Two qubits interact by passing 
through a network of beam splitters and phase shifters along with some ancillary photons. 
The result of this interference will be a very complicated entangled state, but after the 
number of photons in the ancillary modes is measured the state will collapse. Depending 
on the result of the measurement, this process may have performed the desired interaction 
(the controlled sign fiip) or some undesired interaction, in which case the computation has 
failed. Knill, Lafiamme, and Milburn also show that the probability that the measurement 
process reports that the operation has succeeded may be increased arbitrarily close to one by 
increasing the number of ancillary optical modes, photons, beam splitters, etc. This proposal 
keeps many of the advantage of earlier linear optical quantum computers while eliminating 
the exponential growth of the number of optical components with the number of qubits. The 
disadvantage of this scheme is that its calculations succeed only probabilistically, and the 
cost of improving the probability of success is the use of large (but subexponential) numbers 
of ancillary photons and photon detectors. 

The goal of this paper is to investigate the performance of linear optical quantum com- 
puters (LOQC) whose photon detectors have inefficiency and noise. First, we will review 
the LOQC scheme proposed by Knill, Lafiamme, and Milburn. We then discuss a computer 
simulation of LOQC with imperfect detectors. We also examine near deterministic quantum 
teleportation with inefficient detectors. 

I. IDEAL LINEAR OPTICAL QUANTUM COMPUTING 

The goal of LOQC is to build a quantum computer that uses photons to encode informa- 
tion, but does not require any nonlinear medium to allow the photons to interact. Instead, 
we will only use single photon sources, beam splitters, phase shifters, and photon detectors 
which are able to distinguish between 0, 1, 2, ... photons. According to the original scheme 
for LOQC, each qubit is represented by a single photon in a superposition of two optical 
modes. The conventional notation is 

IQ) =a|0i,l2) + 6|li,02). (1) 

Here the logical value of corresponds to the eigenstate that has zero photons in optical 
mode 1 and one photon in optical mode 2, and the logical 1 state has one photon in mode 
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1 and zero photons in mode 2. We will describe a phase shifter using the phase angle 
that it imparts to a single photon, and is the unitary transformation that acts on the 
state of the photons traveling through the phase shifter P^. For example P,/,|0) = |0), and 
P^|l) = e^'^\l). Beam splitters ^e,ii> are similarly described using two angles, 9 and (p. They 
are usually represented by the rotation matrix that transforms the photon creation operators 
ajj into their primed versions in the Heisenberg picture 




cos(^) e'^sm{e) 
-e-^<^sin(^) cos{e) 




(2) 



Because any large unitary matrix can be written as a direct product of two by two unitary 
matrices, we can describe the action of many beam splitters on any number of optical modes 
by giving the matrix that is the direct product of the matrices describing each beam splitter 
and phase shifter. Suppose then that G is such a matrix describing the transformation of 
many optical modes; the elements of G are written as Qij] and G is the unitary operator 
acting on photon states in the Schrodinger picture. To find how a Fock state is transformed 
by G we must compute 



n ^=T {anal + 9i2(A + ■■■+ QimaL) j |0i, O2, 0„) 

i=l y/TT-i- ) 



(3) 



This formalism allows us to generalize the concept of beam splitters to include a much wider 
class of devices, all of which can be constructed by collections of single mode phase shifters 
and conventional two mode beam splitters 0. 

The first element of a LOQC is a device called NS^i that will shift the phase of a mode 
containing two photons, but will leave a single photon undisturbed. NS^i performs the 
transformation 

|Oo) - |Oo) (4) 
|lo)^|lo) (5) 
|2o) ^ -|2o). (6) 

Performing this operation requires the preparation of two ancillary modes containing the 
state |li, O2). The three modes are then sent through a network of beam splitters described 
by the matrix 



( 2-1/4 ^3/V2-2\ 

2-1/4 l/2-l/v/2 



V\/3/v^ 



(7) 



2 I/2-I/V2 V2-I/2 J 



After passing through the network, the state of the three modes will be very complicated, 
but we then measure the number of photons contained in modes 1 and 2. If we measure the 
state 1 11,02), then the operation NS^i was successfully performed. The probability that 
we measure the desired result is 1/4. If some other number of photons is detected then the 
operation has failed, and the information contained in mode has been destroyed. 
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One of the fundamental universal gates for quantum computation is the controlled sign 
flip, which shifts the phase of a two qubit state when both qubits have the logical value of 1 
but does nothing if one of the qubits is 0. In our LOQC notation, this is the transformation 




(8) 
(9) 
(10) 

(11) 



To execute this we must first prepare four ancilla modes containing the state II5, Oe, I7, Os). 
Then we follow the procedure: 

1. Apply the beam splitter to modes 1 and 3. 

2. Perform the NS^i operation on modes 1, 5, and 6. 

3. Perform the NS^i operation on modes 3, 7 and 8. 

4. Apply the beam splitter B_^/4 to modes 1 and 3. 

5. Measure the number of photons in modes 5 through 8. 

If the state II5, Oe, I7, Og) is detected, then the controlled phase shift has been performed 
correctly. If some other state is detected, then the qubits have been destroyed. This will 
succeed with probability (1/4)^ = 1/16, so we call this protocol for the controlled sign flip 
c — Zi/iQ. Soon we will discuss other protocols that will increase the success probability, but 
these are more easily understood in the context of quantum teleportation. 

The goal of quantum teleportation in a LOQC is to transfer the photon state of one 
mode to some other mode without directly connecting the two through a beam splitter. The 
following procedure will transfer the state = a|Oo) + 6|lo) of mode to mode 2. 

1. Prepare the state |li,02). 

2. Apply the beam splitter B^/4_o to modes 1 and 2. This creates the state \ti) — 



3. Apply the beam splitter B_7r/4,o to modes and 1. 

4. Measure the number of photons in modes and 1. 

If the measurement yields |lo, Oi) the mode 2 will contain the state a|0o) — 6|lo), so to return 
this to the original state apply the phase shifter P^r to mode 2. If the measurement 
yields |Oo, li) mode 2 will contain \ip) and nothing must be done. If the measurement has 
some other result, the operation has failed. This procedure will succeed with probability 
1/2, but it can be increased with the following method: 

1. Prepare the state 



l/v^(|0i,l2) + |ll,02)). 



1 



n 



tn) 



+ 1 



^ii)^ior-ioni) 



(12) 



j=0 
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where = i times. |t„) requires 2n optical modes, and they should be 

numbered from left to right in each term of the sum. \tn) can be prepared using only 
beam splitters and the state |1)"'|0)"'. 

2. Operate on modes through n with the beam splitter network Fn, where the elements 
of the matrix describing Fn are {F^jj^k = (e'^'^/*^"+-^^)-'''^/-\/n + 1 

3. Measure the photons in modes through n. 

If the number of photons detected is m, and 1 < m < n, then the teleportation procedure 
has succeeded. The original state occupying mode is found in mode n + m, but this state 
may need to be corrected with a known phase shifter. The probability that 1 < m < n is 
n/ (n + 1). 

To improve the success of the controlled sign flip, we will perform the teleportation 
procedure on each qubit, but rather than use two copies of the state we will use the 
state 

|cs„) = J2 (-l)^""'^^""''^|l)iO)"-'|0)'|l)"-'|l)^'|0)"-^'|0)^'|l)"-^'. (13) 

i,j=0 

Knill, Laflamme, and Milburn give an algorithm for the preparation of this state in [Q. 
According to their algorithm, the preparation of |cs„) requires 6n — 3 controlled phase shift 
operations ||^ each of which succeeds only probabilistically by performing the c — Zi/ie 
operation or by preparing \cSm) for m < n and executing this teleportation method. They 
provide no proof that their method is optimal, and we believe that finding the optimal 
procedure for preparation of |cs„) is a fruitful open problem. Imagine that qubit one occupies 
modes qi and q2, qubit two occupies modes qs and q^, and |cs„) occupies modes 1 through 
An. The algorithm is: 

1. Operate on modes qi and the first n modes of |cs„) with the beam splitter network 
F 

2. Operate on modes q^ and the third n modes of \cSn) with Fn. 

3. Measure modes qi and 1 through n. U 1 < ki < n photons are detected. The contents 
of mode qi are now found in mode n + ki. 

4. Measure modes q^ and 2n + 1 through 3n. If 1 < k2 < n photons are detected, then 
the contents of mode q^ are now found in mode 3n + k2. 

5. Depending on ki and ^2, apply correcting phase shifters to n + ki and 3n + k2. 

Provided that ki and k2 satisfy the above conditions, this procedure will perform the con- 
trolled sign flip on the qubits rather than simply teleporting them. The probability that both 
teleportations succeed is {n/{n + 1))^, and we will refer to this procedure as c — 2;(n/(n+i))2- 
To increase this probability beyond 95% would require n > 39, and to prepare the necessary 
state ICS39) involves 231 uses of the controlled sign flip operation each of which succeeds 
probabilistically. 
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II. SIMULATION OF LOQC WITH IMPERFECT DETECTORS 



In order to investigate the effects of imperfect photon detectors on the operation of the 
LOQC, we have developed a computer simulation using Mathematica. 0] The simulation 
represents the state of the LOQC in the Fock basis in which a basis state assigns an integer 
number of photons to each optical mode. Then superpositions are built out of the basis states 
in the usual manner. The simulation can then calculate how the operation of a collection of 
beam splitters and phase shifters transform the states. To simulate ideal measurement the 
simulation will generate a table displaying the probability to measure all possible numbers of 
photons in each mode being measured, the number of measured photons, and the collapsed 
state produced after the measurement. At this stage the simulation has calculated the 
probabilities that the measured modes of the computer contain 0, 1, 2, or photons, 
we must now calculate the probability that these photons are registered by the detectors. 

We model photon loss using a binomial distribution. |^ If photons enter a detector, 
the probability that n photons are registered by the detector is 



when n < N, a.nd pioss{n, N,l) = when n > N. The inefficiency of the detectors is described 
with the parameter / which is the probability that any given photon will be missed by the 
detector. 

The photon detections may also be corrupted by noise photons. To calculate the proba- 
bility that n noise photons are added to the number of "real" photons each detector registers, 
we apply the Poisson distribution. 



where g is mean number of noise photons each detector is expected to register. Note that 
this Poisson distribution does not represent "quantum noise" in the photon field. Instead, 
we imagine it to represent accidental, for example thermal, excitations inside the detector 
which are independent of the incoming photons. 

The exact physics of the atom-photon interactions in a photon detector are very com- 
plicated, and we make no attempt to model them here. Our goal is only to describe the 
interactions using the two phenomenological parameters I and g. It is difficult to imagine 
how a more detailed understanding of photon detector physics would yield wildly different 
probability distributions, so we are confident that our treatment will give a good under- 
standing of the reliability of LOQC. 

We can apply these probability distributions to the photons entering each detector to 
build a table that shows the probability that each detector registers 1, 2, 3, ... photons and 
the number of "real" photons that arrived at each detector. By adding the probabilities 
for the cases in which the number of photons registered by each detector and the num- 
ber of photons that arrived at each detector are equal to one another and are both equal 
to the number prescribed by the LOQC protocol we can calculate Ps{i,g) the probability 
that a given operation will succeed and the detectors will correctly detect the success. We 
also calculate Pd{i,g), the probability that the detectors report that a given operation has 




(14) 



Pnoise (^) S') 




(15) 
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FIG. 1. (a) shows the probabihty that the NS-i will succeed and the detectors will register 
the success plotted as a function of the detector inefficiency / and the noise level g, for the case 
when the input mode contains two photons, (b) shows the probability that NS-i has failed, but 
the detectors register a success. 

succeeded, when in fact it may have failed, by adding the probabihties that the detectors 
register the number of photons prescribed by the LOQC protocol without regard to the num- 
ber of "real" photons that arrived at the detectors. Using ps and pd, we can also calculate 
Pf{l,g) = 1 — ^ the probability that the detectors register a false positive result, when the 
detectors register the number of photons prescribed by the LOQC protocol but the number 
of photons that arrived at the detectors is not equal to the prescribed number. This is surely 
the most damaging sort of error, because it is disguised as a successful calculation. 

We first applied this analysis to the NS-i operation, which succeeds when one photon is 
detected in mode 1 and zero photons are detected in mode 2. For an ideal device this should 
happen with probability ^. The results of our simulation appear in Fig. ^ where we have 
plotted the ps{l,g) and pf{l,g) for the case when the input mode contains two photons. 

We have also applied a similar simulation to the conditional sign flip c — Zi/ig, when both 
entering qubits are in the logical |1) state. Plots showing the success probability and false 
success probability are shown in Fig. We find that for detectors with / = 0.1 and g = 0.1, 
the operation will succeed in only 3.5% of its operations, and 35% of its apparent successes 
are in fact failures. The probability of obtaining false positive results will be below 1% for 
detectors better than those for which I = g = 0.0025; / = and g = 0.0033; or / = 0.0099 
and g = 0. 

The third system to which we apply our simulation is the controlled sign flip, c — Zi/4, 
in which the state |csi) = I (-|11) + |10) + |01) + |00)) is prepared using a single operation 
of the c — Zi/iQ gate. The two qubits, both in the logical |1) state, are then teleported using 
the |csi) state. Using perfect components, the probability of preparing |csi) is 1/16, and the 
probability of teleporting both qubits is 1/4. Given detectors with I = g = 0.1 we find that 
the probabihty to prepare |csi) is 0.035, slightly more than half of the ideal probability, and 
42% of the apparent successes are failures. Using these detectors the probability to both 
prepare |csi) and to teleport the qubits is 0.0048, only one third of the ideal success rate. 
False positive results can be expected in 58% of the attempts. False positive rates lower 
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FIG. 2. (a) shows the probabihty that the c — -Zi/ig will succeed and the detectors will register 
the success plotted as a function of the detector inefficiency / and the noise level g, for the case 
when both input qubits are in the |1) state, (b) shows the probability that c — -Si/ie has failed, but 
the detectors register a success. 



than 1% can be achieved with detectors better than I = g = 0.0012, / = 0.0020 and g = 0, 
or / = and g = 0.0033. The less ambitious goal of lowering the false positive rate to 10% 
would be possible with detectors better than I = g = 0.013, / = 0.021 and = 0, or / = 
and g = 0.036. These results are plotted in Fig. |^ 



III. NEAR DETERMINISTIC TELEPORTATION WITH IMPERFECT 

DETECTORS 



Because quantum teleportation is such an integral part of the improved versions of the 
controlled sign flip c — Z(^n/{n+i))^ , we sought to understand the rate at which detector related 
errors grows as n the size of the prepared teleportation state increases. Provided that the 
state \tn) has been successfully prepared, we want to examine the probability that the state 
|lo) can be teleported when using detectors with inefficiency /. For simplicity we assume 
that the detectors are noiseless. By examining several cases we have found a general formula 
for Ps{n, I) the probability that the teleportation succeeds and is correctly registered by the 
detectors: 



The probability that exactly i = 1,2,..., or n photons arrive at the detectors is and 
the probability that all of the i arriving photons are registered by the detectors is (1 — /)*. 
According to this expression, for a given detector inefficiency /, there is a critical n = ric for 
which ps is maximized, and increasing n past ric will actually decrease the probability that 
the teleportation will succeed and be correctly registered by the detectors. In Fig. ^ we 
plot the maximum achievable Ps and the ric yielding this probability for a range of detector 
inefficiencies. 
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FIG. 3. (a) shows the probabihty that the c — Zi/4 wiU succeed and the detectors will register 
the success plotted as a function of the detector inefficiency / and the noise level g, for the case 
when both input qubits are in the |1) state, (b) shows the probability that c — has failed, but 
the detectors register a success. 



Of course, any increase in the number of detectors involved in the teleportation effort 
will increase the probability pf{n,l) that the detectors report that the teleportation has 
succeeded, when it has not. To find pf{n,l) we first calculate pd{n,l) the probability that 
the detectors report that the teleportation was successful. This includes all of the events 
counted in Ps{n, I) and all of the ways that photons may be lost so that the detectors register 
1, 2, 3, or n photons. 



Pf{n, I) is then given by 



p,(„.0 = l-^. (18) 

Vdin.l) 



Any increase in n will always increase the probability of receiving false positive results from 
the detectors. In Fig. |^c we have plotted pf{nc,l). To build a teleportation device with a 
success probability of 0.9 requires detectors with / = 0.0055 and n = 19, and for 0.99 success 
rate requires / = 0.00006 and n = 182. Fig. |^ shows how quickly pj increases with n for 
detector inefficiencies of 1% and 0.01%. 

Construction of the near deterministic controlled sign flip requires the preparation of 
the state |cs„), which contains 2n photons distributed throughout An modes. The first 2n 
modes are used to teleport the first qubit and the second 2n modes to teleport the second 
qubit, and after the teleportation the qubits will have received a controlled sign flip. The 
probability that each of these teleportations succeeds is equal to Ps{n, I), and the probability 
that detectors report that each teleportation succeeds is just pd{n,l). The probabilities of 
detecting the various possible numbers of photons during the teleportation of each qubit are 
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FIG. 4. (a) shows the maximum probabihty that the teleportation will succeed and the de- 
tectors will register the success plotted as a function of the detector inefficiency /, for the case 
when the mode to be teleported contains one photon, (b) shows the value of iic that gives the 
maximum success probability, (c) shows the expected rate of false positive results when attempting 
a teleportation using a prepared state with size ric- 
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FIG. 5. This shows the expected rate of false positive results as a function of the size of the 
prepared teleportation state for a device using detectors with inefficiency of 1% in (a) and 0.01% 
in (b). The horizontal axes terminate at ric, the size of the prepared teleportation state giving the 
maximum success probability for each detector. 



independent of one another, so the probabihty to successfully implement the controlled sign 
is {ps{n,l))'^, and the probability of obtaining false positive results is 

For detectors with / = 0.1 the maximum achievable success rate is 0.38, which is obtained 
when n = 4, but 48% of apparent successes are false. To build a controlled sign teleportation 
that succeeds 90% of attempts requires detectors with loss / < 0.0014 and n > 37. A 99% 
reliable gate needs I < 0.000013 and n > 392. 

All of the results in this section rely on the assumption that the detectors produce zero 
noise. We chose to omit noise because of the complexity it adds to the calculations, because 
based on our simulations the errors caused by noise are very similar to those caused by 
loss when both the noise and loss are small, and because of the difficulty of displaying the 
results when varying inefficiency, n, and noise. Knill, Lafiamme, and Milburn mention in 
that one possible way to reduce the occurrence of false positive results in the teleportation 
procedure is to, after measuring k photons in modes through n, also measure the photons 
in modes n+1 through n + k — 1 and modes n + k + 1 through 2n (all of the ancillary modes 
that do not contain the teleported qubit). Then the total number of detected photons should 
equal the number of photons in otherwise one or more has been lost. Although this 
would reduce false positive results it will also reduce the device's success probability because 
sometimes the n — 1 detectors in the second stage will lose photons when the n + 1 primary 
detectors monitoring modes through n function correctly, in which case good events would 
be discarded. If we naively imagine that each detector has the same probability p to lose at 
least one photon, then using this method would reduce the success rate of the teleportation 
from (1 — to (1 — p)^". From all successful teleportations, a fraction of 1 — (1 — p)"^^ 
would be rejected because some of the secondary detectors may malfunction. 
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IV. CONCLUSIONS 



In 1^ Knill, Laflamme, and Milburn conservatively estimate that an accuracy threshold 
better than 99% for controlled sign flip operations is required to construct a reliable LOQC. 
The main result of this paper is the calculation that to achieve this threshold when per- 
forming the c — 2(„/(„+i))2 with a perfectly prepared |cs„) state requires photon detectors 
with inefficiencies lower than / = 0.000013, which are much superior to today's best photon 
detectors. This leads us to search for LOQC schemes that do not rely so heavily on the 
detection of single photons. One such possibility may be the use of coherent state qubits as 
outlined in |TD[ . In their scheme all detectors except one receive a large number of photons, 



and success or failure decisions are based on identifying the detector that registered zero 
photons compared to those detectors that registered a number of photons proportional to 
the (assumed to be large) amplitude of the coherent state. 

This work has made no attempt to analyze the effect of unreliable photon sources and 
beam splitters, which would greatly complicate the analysis. Another open problem is 
the analysis of the many error correction procedures outlined in to understand their 

vulnerability to photon loss and noise. 

Lastly, we should note that this work should not be interpreted as a criticism of the 
many current experimental efforts ITT] to build some basic LOQC elements. Most of these 
experiments apparently plan to use coincidence measurements, which are largely impervious 
to corruption by photon loss or noise because those events can be discarded. However this is 
not a practical strategy to incorporate in a large scale quantum computer, because the qubits 
cannot be measured until the entire calculation involving very large numbers of controlled 
sign flip operations is complete. 
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